The Feigenbaum's 5 for a high dissipative bouncing ball model 
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We have studied a dissipative version of a one-dimensional Fermi accelerator model. The dynamics 
of the model is described in terms of a two-dimensional, nonlinear area-contracting map. The 
dissipation is introduced via innelastic collisions of the particle with the walls and we consider the 
dynamics in the regime of high dissipation. For such a regime, the model exhibits a route to chaos 
known as period doubling and we obtain a constant along the bifurcations so called the Feigenbaum's 
number S. 
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I. INTRODUCTION 

The origin of cosmic radiation and a mechanism 
trough which it could acquire enormous energies has in- 
trigued many physicists and mathematicians along the 
last decades. In particular, in the year of 1949 En- 
rico Fermi [ij proposed a very simple model that qual- 
itatively describes a process in which charged particles 
were bounced via interaction with moving magnetic fields. 
This heuristic idea was then later modified to encompass 
in a suitable model that could give quantitative results 
on the original Fermi's model. After that, many different 
versions of the model were proposed and studied in 
different context and considering many approaches and 
modifications One of them is the well known 

Fermi- Ulam model (FUM). Such model consists of a clas- 
sical point-like particle moving between two walls in the 
total absence of any external field. One of the walls is 
considered to be fixed while the other one moves period- 
ically in time. Despite this simplicity, the structure of 
the phase space is rather complex and it includes a large 
chaotic sea surrounding KAM islands and a set of invari- 
ant spanning curves limiting the energy of a bouncing 
particle. However, the introduction of innelastic colli- 
sions on this model is enough to destroy such a mixed 
structure and the system exhibits attractors. Depend- 
ing on the initial conditions and control parameters, one 
can observe a chaotic attractor characterised by a posi- 
tive Lyapunov exponent. By a suitable control parameter 
variation, the chaotic attractor might be destroyed via a 
crisis event p"3 - |20l |. After the destruction, the chaotic 
attractor is replaced by a chaotic transient [l9j . 

In this paper, we revisit a dissipative version of the 
Fermi-Ulam model seeking to understand and describe 
the behaviour of the dynamics for the regime of high 
dissipation. Moreover, we investigate a phenomenon 
quite common observed for a large variety of dissipative 
and nonlinear systems known as period doubling route to 
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chaos. This route shows a sequence of doubling bifur- 
cations connecting regular-periodic to chaotic behaviour. 
The main goal of this paper is then to examine and show 
that the Fermi-Ulam model obey the same convergence 
ratio as that obtained by Feigenbaum jUj for the dynam- 
ics of the logistic map. Such a constant value has been 
further called as "The Feigenbaum's S" . 

The organisation of this paper is as follows: In sec- 
tion |lT] we present all the details needed for the construc- 
tion of the nonlinear mapping. Our numerical results 
are discussed in section IIIII and our conclusions and final 
remarks are drawn in section ITVl 



II. A DISSIPATIVE FERMI-ULAM MODEL 

Let us describe the model and obtain the equations of 
the mapping. The model we are dealing with consists of a 
classical point-like particle confined in and bouncing be- 
twen two rigid walls. One of them is assumed to be fixed 
at the position x = I and the other one moves periodically 
in time according to x^it) = ecos(wt). Thus, the mov- 
ing wall velocity is given by Vwit) = —eujshi{ujt) where e 
denotes the amplitude of oscillation and lo is the angular 
frequency of the moving wall, respectively. Aditionally, 
the motion of the particle does not suffer infiuence of 
any external field. We assume that all the collisions with 
both walls are inelastic. Thus, we introduce a restitu- 
tion coefficient for the fixed wall as a G [0, 1] while for 
the moving wall we consider /3 e [0, 1]. The completely 
inelastic collision happens when a = /3 = and that a 
single collision is enough to terminate all the particle's 
dynamics. On the other hand, when a = (3 = 1, corre- 
sponding to a complete elastic coUision^ all the results for 
the nondissipative case are recovered [17]. 

We describe the dynamics of the particle using a map 
T where the dynamical variables are (vmtn) and Vn is 
the velocity of the particle at the instant i„. The in- 
dex n denotes the n*^ collision of the particle with the 
moving wall. Starting then with an initial condition 
(vmtn), with initial position of the particles given by 
Xp{t) = ecos(wt„) with Vn > 0, the dynamics is evolved 
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and the map T gives a new pair of {vn+i,tn+i) at the 
(n + 1)*'' collision. It is important to say that we have 
three control parameters, namely: e, / and w and that the 
dynamics does not depend on all of them. Thus it is con- 
venient to define dimensionless and more appropriated 
variables. Therefore, we define e = e/Z, Vn = Vn/ujl and 
measure the time in terms of the number of oscillations 
of the moving wall, consequently 0„ = Lotn. Incorporat- 
ing this set of new variables into the model, the map T 
is written as 



T 



K+i =Kr-(l + /3)esin(0„+i) 
4>n+i ^ 4>n + Ar„ mod(27r) 



(1) 



where the expressions for both V* and AT„ depend on 
what kind of collision occurs. There are two different 
possible situations, namely: (i) multiple collisions with 
the moving wall and, (ii) single collision with the moving 
wall. Considering the first case, the expressions are V* = 
—PVn and Ar„ = (pc, where 0c is obtained as the smallest 
solution of G{(t>c) — with (f)c G (0,27r]. A solution for 
G{(j)c) = is equivalent to have that the position of the 
particle is the same as that of the moving wall at the 
instant of the impact. The function G(0c) is given by 



G(0c) 



-I- 0c) - ecos((/)„) - V„(j)c 



(2) 



If the function G{(pc) does not have a root in the interval 
(/)c G (0, 2tt], we can conclude that the particle leaves the 
collision zone (the collision zone is defined as the interval 
X G [— e, e]) without suffering a successive collision. 

Considering now the case (ii), i.e. the case where the 
particle leaves the collision zone, the corresponding ex- 
pressions are V* — —a(3Vn and Ar„ — + (pr + (ph with 
the auxiliary terms given by 



1 - £COs(0n) 
Vn 



1 - e 
aVn 



(3) 



Finally, the term (pc is obtained as the smallest solution 
of F{(j)c) = for 0c G [0, 27r). The expression of -F(0c) is 
given by 



+ (j)r + 4>l + (fic) ~ £ + aVn4>c 



(4) 



Both, Eq. ([2]) and Eq. Q must be solved numerically. 

After some algebra it is easy to shown that the deter- 
minant of the Jacobian Matrix for the case (i) is given 

by 



detiJ) = 13^ 



Vn+i + esin(0„+i) 
while for the case (ii) we obtain 

Vn + esin(0„) 



det{J) = a^l3^ 



Vn^ 



e sm 



>n+l 



(5) 



(6) 



We therefore can conclude based on the above result that 
area preservation is obtained only for the case oia — (3 — 
1. 



III. NUMERICAL RESULTS 

In this section we discuss some numerical results con- 
sidering the case of high dissipation. We call as high 
dissipation the situation in which a particle loses more 
than 50% of its energy upon colision with the moving 
wall. Moreover, we have considered as fixed the values 
of the damping coefficients a = 0.85 and /? = 0.50. To 
illustrate that the dynamics of the system exhibits dou- 
bling bifurcation cascade, it is shown in Fig. [TJa) the 
behaviour of the asymptotic velocity plotted against the 
control parameter e, where the sequence of bifurcations 
is evident. A similar sequence is also observed for the 
asymptotic variable 0, as it is shown in Fig. [IJb). Note 
that all the bifurcations of same period in (a) and (b) 
happen for the same value of the parameter e. At the 
point of bifurcations, one can also observe that the posi- 
tive Lyapunov exponent shows null value, as it is shown 
in Fig. [Ijc). It is well known that the Lyapunov ex- 
ponents are important tool that can be used to classify 
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FIG. 1: Bifurcation cascade for (a) 4> and (b) V both plotted 
against e. In (c) it is shown the Lyapunov exponent associ- 
ated to (a) and (b). The damping coefficients used for the 
construction of the figures (a), (b) and (c) were a = 0.85 and 
/3 = 0.50. 
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orbits as chaotic. As discussed in ^2^, the Lyapunov 
exponents are defined as 

A, = hm - In I A, I , j = l,2 , (7) 

n— ^oo fi 

where Aj are the eigenvalues of M — n"=i ^ii^ij 4'i) ^^id 
Ji is the Jacobian matrix evaluated over the orbit {Vi, 4'i). 
However, a direct implementation of a computational al- 
gorithm to evaluate Eq. ([7]) has a severe limitation to 
obtain M . Even in the limit of short n, the compo- 
nents of M can assume very different orders of magni- 
tude for chaotic orbits and periodic attractors, yielding 
impracticable the implementation of the algorithm. In 
order to avoid such problem, we note that J can be writ- 
ten as J = QT where is an ortoghonal matrix and 
T is a right up triangular matrix. Thus we rewrite M 
as M = J„ J„_i . . . J2QiQi^Ji, where Ti = Q^^Ji- A 
product of defines a new Ji^. In a next step, it is 

easy to show that M = JnJn-i ■ ■ ■ J3Q2'd2^ -^2^1- The 
same procedure can be used to obtain T2 = J2 
so on. Using this procedure, the problem is reduced to 
evaluate the diagonal elements of Ti : Tii,Tj2. Finally, 
the Lyapunov exponents are given by 

n 

X = hm -V ln|r;,| , j = 1,2 . (8) 

If at least one of the Aj is positive then the orbit is 
classified as chaotic. We can see in Fig. ([T]) (c) the be- 
havior of the Lyapunov exponents corresponding to both 
Figs. ([T])(a,b). It is also easy to see that when the bifur- 
cations happen, the exponent A assumes the zero value 
at same values of the control parameter e where the bi- 
furcation hold. Another interesting observation is that 
the Lyapunov exponent, in some regions, assumes a con- 
stant and negative values, such behavior occurs because 
the eigenvalues of the Jacobian Matrix become complex 
numbers (imaginary numbers)). 

Let us now use the points where the Lyapunov expo- 
nent assumes the null value to characterise a very inter- 
esting property of the model. As it was shown by Feigen- 
baum [2l|, along the bifurcations the dynamics exhibits 
an universal feature. It implies that all the bifurcations 



happen at the same rate of convergence for the bifurca- 
tion diagram changing from periodic to chaotic behavior. 
The procedure used to obtain the Feigenbaum value 5 
consists of: let ei represents the control parameter value 
at which period-1 gives birth to a pcriod-2 orbit, 62 is the 
value where period-2 changes to period-4 and so on. In 
general the parameter e„ corresponds to the control pa- 
rameter value at which a period-2" orbit is born. Thus, 
we write the Feigenbaum's 5 as 

5„ = hm (9) 

The numerical value for the constant number 5 ob- 
tained by Eq. ^ and considering a sufficient large values 
for n is 5 = 4.66920161 . . .. Considering the numerical 
data obtained through the Lyapunov exponents calcula- 
tion, the Feigenbaum's 5 obtained for the Fermi-Ulam 
model is (5 = 4.64(4). We have considered in our sim- 
ulations only the bifurcations of fourth to eighth order 
since the numerical results are very hard to be obtained 
for higher orders in the bifurcations. Our result is in a 
good agreement with the Feigenbaum's universal 5. 



IV. CONCLUSION 

As a summary, we have studied a dissipative version 
of the Fermi-Ulam model. We introduce dissipation into 
the model through inelastic collisions with both walls. 
We have shown that for regimes of high dissipation, the 
model exhibits a sequence of doubling bifurcation cas- 
cade. For this cascade, we obtained the so called Feigen- 
baum's number &s 5 — 4.64(4). 
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